常微分方程技术及其在固体力学计算中的应用

常微分方程(ODE——Ordinary Differential Equation)边值问题的最新计算求解技术的迅速发展推出了一批高质高效的通用软件,而工程中大量的ODE问题并非呈现为这些求解器(Solver)所接受的标准形式。然而,运用一些简单的ODE变换技巧可以将大量的不同类型的特殊问题转化为标准形式。本文列举了若干常用的变换技巧,并广泛地应用于各种固体力学问题的计算中,使大量的ODE问题在形式上得到统一,得以用标准的ODE Solver方便有效地求解。     

Bifurcations et solutions multiples en cavité 3D différentiellement chaufféeBifurcations and multiple solutions in a differentially heated cubic cavity

In this paper, a differentially heated square/cubic cavity is studied by performing three-dimensional direct numerical simulations. The first bifurcation observed at $\text{R}{}_{\mathrm{a}}\text{≈3.2×10}{}^7$ is due to the 3D vortex structures generated at the end regions of vertical boundary layers near the median plane. The main results of this Note are that the flow returns to a steady state for higher values of the Rayleigh number $\text{R}{}_{\mathrm{a}}$ (7×10⁷ and 10⁸ for example) still exhibiting these 3D vortex structures, and that multiple steady flows which differ by their symmetry properties, are obtained for $\text{R}{}_{\mathrm{a}}\text{=10}{}^8$. However, the flow reverts to unsteadiness for $\text{R}{}_{\mathrm{a}}\text{=3×10}{}^8$. In this latter case, the instability is due to the vertical boundary layers. To cite this article: G. de Gassowski et al., C. R. Mecanique 331 (2003).     

On the fundamental sloshing frequency in the ice-fishing problemSur la fréquence de ballottement principale dans le problème de la glace avec trous pour pêcher

The sloshing problem is considered in a half-space covered by a rigid dock with apertures. The dependence of the fundamental sloshing frequency on the shape of the free surface region is studied. It is proved that the inequality holds between the fundamental eigenvalues corresponding to two different regions if some conditions are fulfilled. These conditions are verified for particular classes of regions of a fixed area in order to demonstrate that the disk yields the maximum of the fundamental eigenvalue for each of these classes. On the other hand, examples of regions are constructed for which the fundamental eigenfrequency is larger than that for the circular aperture of the same area and even as large as one wishes. To cite this article: V. Kozlov, N. Kuznetsov, C. R. Mecanique 330 (2002) 723–728.     

On the semiclassical density of states of quasi-integrable mechanical systems

For a class of quasi-integrable mechanical systems, we show that in the semiclassical limit the Bohr-Sommerfeld quantization applied to successive truncations of the Birkhoff series and a power series expansion of Weyl's formula give the same asymptotic number of states below a given energy.

---

Sommario Si mostra che per una classe di sistemi quasi integrabili, nel limite semiclassico, la quantizzazione di Bohr-Sommerfeld applicata a successivi troncamenti della serie di Birkhoff e l'espansione in serie di potenze della formula di Weyl danno lo stesso andamento asintotico per la densita' degli stati integrata.

     

A centrally cracked thin circular disk, Part I: 3D Elastic–plastic finite element analysis

A full field solution, based on small deformation, three-dimensional elastic–plastic finite element analysis of the centrally cracked thin disk under mode I loading has been performed. The solution for the stresses under small-scale yielding and lo!cally fully plastic state has been compared with the HRR plane stress solution. At the outside of the 3D zone, within a distance of rσ_o/J=18, HRR dominance is maintained in the presence of a significant amount of compressive stress along the crack flanks. Ahead of this region, the HRR field overestimate the stresses. These results demonstrate a completely reversed state of stress in the near crack front compared to that in the plane strain case. The combined effect of geometry and finite thickness of the specimen on elastic–plastic crack tip stress field has been explored. To the best of our knowledge, such an attempt in the published literature has not been made yet. For the qualitative assessment of the results some of the field parameters have been compared to the available experimental results of K, gives a fair estimate of the crack opening stress near the crack front at a distance of order 10⁻² in. On the basis of this analysis, the Linear Elastic Fracture Mechanics approach has been adopted in analyzing the fatigue crack extension experiments performed in the disk (Part II).     

The effect of velocity on rigid wheel performance

A simulation model to predict the effect of velocity on rigid-wheel performance for off-road terrain was examined. The soil–wheel simulation model is based on determining the forces acting on a wheel in steady state conditions. The stress distribution at the interface was analyzed from the instantaneous equilibrium between wheel and soil elements. The soil was presented by its reaction to penetration and shear. The simulation model describes the effect of wheel velocity on the soil–wheel interaction performances such as: wheel sinkage, wheel slip, net tractive ratio, gross traction ratio, tractive efficiency and motion resistance ratio. Simulation results from several soil-wheel configurations corroborate that the effect of velocity should be considered. It was found that wheel performance such as net tractive ratio and tractive efficiency, increases with increasing velocity. Both, relative wheel sinkage and relative free rolling wheel force ratio, decrease as velocity increases. The suggested model improves the performance prediction of off-road operating vehicles and can be used for applications such as controlling and improving off-road vehicle performance.     

岩石（煤）的碎胀与压实特性研究

本文通过自行设计的一套岩石（煤）的松散和压实试验方案及相应的实验设备，对兖州矿区岩（煤）样进行了较为系统的碎胀与压实特性研究，测定了岩石（煤）的碎胀系数、碎胀曲线、压实曲线和侧压曲线．     

A highly redundant coordinate formulation for constrained rigid bodies

The paper presents a formulation for the dynamic analysis of rigid multibodies. An introductory part carries out the kinematic analysis and the definition of the highly redundant differential framework along with the choice of unknowns and equations. From the differential formulation the variational principles, either in Lagrangian or Hamiltonian form, are developed. The Hamiltonian formulation is then used to develop the numerical approximation by applying the finite element method in time. The application of the method in its multifield form is discussed and a solution algorithm is proposed. Some examples are finally presented in order to verify the effectiveness of the formulation.

---

Sommario Il lavoro presenta una formulazione per lo studio della dinamica dei sistemi multicorpo rigidi. Nella parte introduttiva viene svolta l'analisi cinematica e si definisce il quadro differenziale con la scelta delle incognite e delle equazioni. Dalla formulazione differenziale vengono poi sviluppati dei principi variazionali nella forma Lagrangiana ed Hamiltoniana. La formulazione Hamiltoniana é quindi utilizzata per sviluppare l'approssimazione numerica col metodo degli elementi finiti di tempo. Viene discussa l'applicazione del metodo nella forma multi-campo e viene proposto un algoritmo di soluzione. Da ultimo, vengono discussi alcuni esempi per verificare la correttezza della formulazione.

     

A new two-point deformation tensor and its relation to the classical kinematical framework and the stress concept

Starting from the issue of what is the correct form for a Legendre transformation of the strain energy in terms of Eulerian and two-point tensor variables we introduce a new two-point deformation tensor, namely H=(F−F^−T)/2, as a possible deformation measure involving points in two distinct configurations. The Lie derivative of H is work conjugate to the first Piola–Kirchhoff stress tensor P. The deformation measure H leads to straightforward manipulations within a two-point setting such as the derivation of the virtual work equation and its linearization required for finite element implementation. The manipulations are analogous to those used for the Lagrangian and Eulerian frameworks. It is also shown that the Legendre transformation in terms of two-point tensors and spatial tensors require Lie derivatives. As an illustrative example we propose a simple Saint Venant–Kirchhoff type of a strain-energy function in terms of H. The constitutive model leads to physically meaningful results also for the large compressive strain domain, which is not the case for the classical Saint Venant–Kirchhoff material.     

双周期裂纹压电材料反平面剪切问题断裂力学分析

研究压电材料双周期裂纹反平面剪切与平面电场作用的问题．运用复变函数方法，获得了该问题严格的闭合解，并由此给出了裂纹尖端应力强度因子和电位移强度因子的精确公式．数值算例显示了裂纹分布特征对材料断裂行为的重要影响．叠间小裂纹能够对主裂纹的应力和电位移场起着屏蔽作用，相反行间小裂纹却起着放大作用，至于钻石形分布裂纹的影响规律则更为复杂．对于某些特殊情形给予了解答并导出一系列有意义的结果。